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Testing the 'dark-energy '-dominated cosmology via the 
Solar-system experiments 
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ABSTRACT 

The effect of 'dark energy' (i.e. the A-term in Einstein equations) is sought for at 
the interplanetary scales by comparing the rates of secular increase in the lunar or- 
bit obtained by two different ways: (1) measured immediately by the laser ranging 
and (2) estimated independently from the deceleration of the Earth's proper rota- 
tion. The first quantity involves both the well-known effect of geophysical tides and 
the Kottler effect of A-term (i.e. a kind of the 'local' Hubble expansion), while the 
second quantity is associated only with the tidal influence. The difference between 
them, 2.2±0.3 cmyr , can be attributed just to the local Hubble expansion with 

rate fl£ oc) = 56±8 kms" 1 Mpc" 1 . Ass uming that Hubble expansion is formed lo- 
cally only by the uniformly distributed dark energy (A-term), while globally also by 
a clumped substance (for the most part, the cold dark matter), the total (large-scale) 
Hubble constant should be Hq = 65±9 kms" 1 Mpc" 1 . This is in reasonable agreement 
both with the commonly-accepted WMAP result, Hq = 71±3.5 kms -1 Mpc -1 , and 
with the data on supernovae la distribution. The above coincidence can serve as one 
more argument in favor of the dark energy. 
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1 INTRODUCTION 

According to the recent astronomical data, the most part of 
energy density in the Universe (up to 75%) is in the 'dark' 
form (such as the so-called 'quintessence', inflaton poten- 
tial, polarization of vacuum and so on), which is effectively 
de scribed by A-term in t he Einstein equations (e.g. r eviews 
bv lvan den Berghlll999t IcherninlEoOlt lKrausJl20ol etc.). 
All arguments in favor of the dark energy were obtained so 
far from the observational data related to very large (in- 
tergalactic) scales, such as a distribution of supernovae la 
as function of their redshift, the spectrum of fluctuations of 
the cosmic microwave background radiation, the spectra of 
absorption lines from distant sources, or Lya forest, and so 
on. 

Is it possible to find a manifestation of the dark en- 
ergy at much less scales (e.g. inside the Solar system)? In 
general, such effects can be expected from the solution of 
the equations of General Relativity for a point-like mass M 
in the A- dominat e d (de Sitter) Universe, which was ob- 
tained by iKottlerl 1^918) a very short time after the orig- 
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inal Schwarzschild solution. (More details are given in Ap- 
pendix [X]) 

The presence of A-term should change, particu- 
larly, the standard relativistic shift of Mercury's perihe- 
lio n, predicted by General Relativity. This was the idea 
bv lCardona fc Teieirol 1^)98), who proposed using the mea- 
sure of the uncertainty in our knowledge of Mercury's per- 
ihelion shift to impose the upper bound on A. The result 
obtained was not so good as other cosmological estimates 
but, surprisingly, the accuracy was worse by only 1^-2 or- 
ders of magnitude. So, according to the above-cited authors, 
improvement of the value of the shift by one or two decimal 
digits should make such method of determination of A com- 
petitive with the observations at large scales. A more skep- 
tic al viewpoint on the same subject was presented recently 
bv llorid feOOrjft . 

In any case, since accuracy of the above metho d is still 
insufficient , it was proposed in our previous papers llDuminl 
l200ll 120031) to utilize the data of radial (rather than angu- 
lar) measurements of the Moon to reveal anomalous increase 
in its orbit produced by the A-term in metric HA4I - IA8L 
under assumption that the central mass M belongs to the 
Earth. This looks formally as 'local' Hubble expansion. Un- 
fortunately, the result of the earlier works was quite strange: 
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Figure 1. Schematic illustration of the Einstein— Straus theorem. 



the 'local' Hubble constant was found to be approximately 
in the middle between zero and the standard intergalactic 
value, which did not allow a reasonable quantitative inter- 
pretation. 

The aim of the present work is to describe a much im- 
proved analysis of the available observations and to show 
that its results are in reasonable theoretical agreement with 
the large-scale data. 



2 THEORETICAL VIEWPOINTS ON THE 
PROBLEM OF LOCAL HUBBLE 
EXPANSION 

The secular increase in planetary radii due to the A-term, 
discussed in Appendix^ formally looks like a kind of the 
Hubble expansion. So, let us briefly discuss why is it nec- 
essary to reexamine the problem of local Hubble expansion 
just in the context of 'dark-energy'-dominated cosmological 
models? 

In general, Hubble dynamics at the small scales is 
stu died for a long t ime, starting from the pioneering work 
by iMcVittiel il933f) . Although the results by various au- 
thors ob tained by now were quite contradictory (e.g. re- 
view by iBonnorl 1200(1 and references therein), the most 
popular point of view was that the Hubble expansion man- 
ifests itself only at the sufficiently large distances (from 
a few Mpc) and is absent at t he less scales at all (e.g. 
iMisner. Thorne fc Wheelerlll973h . There were a few argu- 
ments in favor of such conclusion. 

The first of t hem is based on the so- called Einstein- 
Straus theorem feinstein fe Straus! Il945fl . Let us con- 
sider the Friedmann-Robertson-Walker (FRW) cosmolog- 
ical model, uniformly filled with some kind of matter. Next, 
let us assume that the background matter inside a specified 
sphere is cut off and concentrated in the point in its centre, 
as illustrated in Fig. Q Then, using equations of General 
Relativity, it can be shown that gravitational field inside the 
empty cavity is described by the purely static Schwarzschild 
metric and begins to experience the cosmological expansion 
only outside its contact boundary with the background mat- 
ter distribution; this boundary moving just with the Hubble 
velocity. Therefore, there is no any Hubble expansion in the 
vicinity of the mass M. 

Unfortunately, despite a mathematical elegance of this 
result, it is absolutely unclear how can the above-mentioned 
cavities be identified for the real astronomical objects. More- 
over, this theorem becomes evidently inapplicable at all in 
the case of 'dark-energy'-dominated cosmology, because it 



is meaningless to consider an empty cavity in the vacuum 
energy distribution. 

The second argument against the Hubble expansion at 
small scales is based on a quasi-Newtonian treatment of 
Hubble effect in a small volume as a tidal-like action by 
dista nt matter (e.g. the recent work bv lDommguez fc Gaitel 
l200ll and references therein). The final conclusion usually 
derived by this way is that there should be no Hubble ex- 
pansion in the gravitationally-bound systems (i.e. the ones 
whose kinetic energy is less than potential energy), such as 
the planetary systems or stars in galaxies. According to this 
criterion, the Hubble expansion should manifest itself only 
from 5^-10 Mpc, as was actually observed for a long time. 

Unfortunately, the tidal treatment of Hubble effect is 
not well justified from the theoretical point of view. Be- 
sides, more accurate astronomical measurements in the re- 
cent few years revealed a well-formed Hubble flow down to 
the distances of about f Mpc, which are c onsiderably less 
than the scale of gravitational unbinding dCherninl l200ll : 
Ekholmetal ] |200lTl . Moreover, in the A-dominated cosmol- 
ogy the tidal effects cannot be of primary significance just 
because of the perfectly uniform distribution of the dark 
energy and, therefore, may be expected only from minor 
constituents of the Universe. 

At last, one more approach for treating the influence of 
cosmic expansion on the dynamics of small-scale systems is 
based on the Einstein-Infeld-Hoffmann (EIH) surface inte- 
gral method, which enables to derive the equations of motion 
of the particles immediately from the field equations. This 
method is based on the integration of field equations over 
the small closed surfaces surrounding the point-like sources 
of the field and subsequent use of independence of these in- 
tegrals on the particular surfaces. Such approach is widely 
employed to obtain the post-Newtonian equations of mo- 
tion against the background of flat (Minkowski) space-time 
and wa s applied a l so to t he problem of local Hubble expan- 
sion by lAndersonl JT995). It was found that planetary sys- 
tems should really expand but with a rate much less than 
for the entire Universe. Unfortunately, the EIH method also 
becomes inappropriate if the space is filled everywhere with 
the dark energy, because the surface integrals are no longer 
invariant when the integration surfaces are moved. 

Finally, a frequent experimental argument against 
Hubble expansion within the Solar system is based on 
the available constraint on time variation in the grav- 
itational constant derived from the planetary dynamics, 
which is now as strong as G/G = (4±9)xl0~ 13 yr _1 
iWilliams. Turvshev fc Boggsl I2OO4I) . So, as was concluded 
by these authors, "the G/G uncertainty is 83 times smaller 
than the inverse age of the Universe, to = 13.4 Gyr . . . Any 
isotropic expansion of the Earth's orbit which conserves an- 
gular momentum will mimic the effect of G on the Earth's 
semimajor axis, a/a — —G/G . . . There is no evidence for 
such local (~1 AU) scale expansion of the solar system." 

Unfortunately, the above-stated equivalence between 
the effect of variable G and the cosmological expansion is 
based solely on the Newtonian arguments. A more accurate 
treatment of this problem in the framework of General Rel- 
ativity for a general case of the multi-component Universe 
is very difficult. Nevertheless, it can be performed in ana- 
lytic form for the particular case of a point-like mass in the 
Universe filled only with dark energy (the A-term), and the 
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Table 1. Rates of secular increase in the mean Earth-Moon distance. 



Method 


Immediate measurement by 


Independent estimate from the 




the lunar laser ranging 0- 


Earth's tidal deceleration'' 


Effects involved 


(1) geophysical tides 


(1) geophysical tides 




(2) local Hubble expansion 




Numerical value 


3.8±0.1 cmyr" 1 


1.6±0.2 cmyr" 1 



a The data were taken from the review bv lDickev et alj Jl994l) . 

^ Formula Q was used with the rate of Earth's diurnal deceleration (S) . 



corresponding results are outlined in Appendix [X] As fol- 
lows from expressions (IA9II - 1A12L the A-dependence of a 
few terms of the resulting metric tensor can really be rein- 
terpreted as the effect of variable G, but this is not true in 
general: there are some terms whose A-dependence is radi- 
cally different from any variations in G. 



3 ANALYSIS OF OBSERVATIONAL DATA 

Since all the commonly-used arguments against the small- 
scale Hubble expansion fail in the case of dark energy, it be- 
comes reasonable to seek for the corresponding effect; and 
the most sensi tive tool seems to be the lunar las er rang- 
ing (LLR) (e.g. lDickev et al.lll994UNordtvedtll999l) . For ex- 
ample, if we assume that planetary systems experience the 
Hubble expansion with the same rate as everywhere in the 
Universe (60-^70 kms" 1 Mpc" 1 ), then average radius of the 
lunar orbit R should increase by approximately 50 cm for the 
period of 20 years. On the other hand, the accuracy of LLR 
measurements durin g the last 20 years was maintained at 
the level of 2-^3 cm JPickev et ai]ll994) : so the perspective 
of revealing the local Hubble effect looks very good. 

The main obstacle that needs to be got around is to ex- 
clude the effect of geophysical tides, which also contributes 
to the secul ar increase in the Earth-Moon distance. As is 
known (e.g. lKaulalfl968l) . because of dissipative effects, the 
Earth's tidal bulge, formed by the lunar attraction, is not 
perfectly aligned in the direction to the Moon but slightly 
shifted towards the Earth's proper rotation. Therefore, there 
is a torque moment, which decreases the proper angular mo- 
mentum of the Earth and increases the orbital momentum of 
the Moon. As a result, the average Earth-Moon distance R 
gradually increases in the course of time. 

According to the law of conservation of angular momen- 
tum, the secular variation in R is related to the change in 
the Earth's diurnal period Te by the simple formula: 

R = kfv, (1) 

where fe = 1.81xl0 5 cms -1 (a more detailed discussion ca n 
be found, for example, in our previous work. lDumm1 E003'l. 
So, if Te is known from independent astrometric measure- 
ments of the Earth's rotation deceleration with respect to 
distant objects, then relation Q can be used to exclude 
the influence of geophysical tides and, thereby, to reveal a 
probable presence of local Hubble expansion. 

Unfortunately, Te is not a well-defined quantity: its val- 
ues derived from the observations in telescopic era and from 



the various sets of pre-telescopic data appreciab ly differ from 
each other (e.g. IStephenson fc Morrisor]|l984h . It is not so 
important for us now what is the reason for these discrep- 
ancies: this may be the systematic errors in the earliest data 
or, for example, the tectonic processes appreciably changed 
the Earth's moment of inertia just before the period of tele- 
scopic observations. Since LLR data refer to the last decades, 
they should be confronted with the most recent values of the 
Earth's rotation deceleration. The corresponding telescopic 
data, starting from the middle of the 17th century, were 
processed by a few researches; and one of the most detaile d 
compilations was presented recently bv ISidorenkovl J2002). 

Of course, the value of secular trend derived from the 
quite short time series can suffer from the considerable peri- 
odic and quasi-periodic variations in Te. So, the main aim of 
our statistical analysis outlined in Appendix iBl was to esti- 
mate as carefully as possible the 'mimic' effect of such vari- 
ations. Taking into account the corresponding uncertainty, 
the resulting value can be written as 

T E = (8.77±1.04)xl0~ 6 syr" 1 . (2) 

The results of the entire analysis of LLR vs. the astro- 
metric data are summarized in Tabled The excessive rate of 
increase of the lunar orbit 2.2±0.3 cmyr -1 can be attributed 
to the local Hubble expansion with rate 

H^ loc) = 56±8 kms" 1 Mpc" 1 . (3) 

This value is appr eciably greater than in our earlier 
work jDuminl l2003tl . where it was found to be only 
33±5 kms" 1 Mpc" 1 . This was because of using a sub- 
stantially different rate of the Earth's deceleration, Te = 
1.4xl0~ 5 syr" 1 . The last-mentioned valu e was o btained for 
the first time bv IStephenson fc Morrison! jf?)* 1 i. who used 
the telescopic observations supplemented by a much longer 
series of medieval Arabian data; and their result was inac- 
curately cite d in a number of subsequent reviews and mono- 
graphs (e.g. |Pertsevll2000t) as derived from the telescopic 
observations alone. 

Finally, it should be mentioned that the basic rela- 
tion was written under assumption of constant moment 
of inertia of the Earth, wh ich is a ques t ionab le item. For 
example, over 20 years ago lYoder et alj (I1983T) found that 
Earth's oblateness, commonly characterized by the gravi- 
tational harmonic coefficient J2, was decreasing. This was 
interpreted as viscous rebound of the solid Earth from the 
decrease in load due to the last deglaciation. The observed 
secular effect J2 = — 3xl0~ n yr" 1 resulted in the Earth's 
spin-up due to decreasing moment of inertia and, thereby, 
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enabled the above-cited authors to get a reasonable agree- 
ment between the various sets of data on the Earth's rota- 
tion. Therefore, following this approach, it would be unnec- 
essary to take into consideration any other influences, such 
as the cosmological Hubble expansion. 

Unfortunately, the early results bv lYoder et al.l <ll983t) 
were not confirmed by the most recent studies. For ex ample, 
as follows from the analysis bv lBourda fc Capitaind (I2004T) 
performed over a sufficiently long time interval (1985-2002), 
the coefficient J2 has a much smaller secular trend but a 
considerable oscillatory component with a period of two 
decades. An even more striking disagreeme nt with the earlier 
data was obtained bv lCox fc ChaoN2002T) . who found that 
since 1997 or 1998 the secular trend in J2 has approximately 
the same absolute value as reported bv lYoder et al.l dl983l) 
but the opposite sign (namely, J2 = +2.2xl0 -11 yr _1 ). 
Therefore, we can conclude that (1) a considerable disagree- 
ment between the LLR and astrometric data still exists and 
(2) the coefficient J2 experiences most probably a quasi- 
periodic variation with a typical time scale of a few decades. 
The last-mentioned property justifies usage of formula Q, 
because the temporal variations in the Earth's moment of 
inertia should be averaged out when data on the Earth's 
rotation are taken for the period of ~350 years. 



4 THEORETICAL INTERPRETATION 

How can the value © be interpreted? ft is reasonable 
to assume that local Hubble expansion is formed only by 
the uniformly-distributed dark energy (A-term), while the 
irregularly-distributed (aggregated) forms of matter begin 
to affect the Hubble flow at the larger distances, thereby 
increasing its rate up to the standard intergalactic value. 

If the universe is spatially flat and filled only with vac- 
uum and a dust-like ('cold') matter, with densities pao and 
Pdo respectively, then 



Ho 



8ttG 



V PAO + pDO 



(4) 



(e.g. lLandau fc Lifshit dll975l) . So, if Ho is formed locally 
only by pao, while globally by both these terms, pao and 
Pdo (or, in terms of the relative densities, Qao = Pao/ Per 
and £2do = Pdo/Pci), then 



H, 



(loc) 



Ho 



1 + 



-1/2 



(5) 



At the commonly-accepted values Oao = 0.75 and £2do : 
0.25, we get H /H { loc) ^ 1.15. Therefore, 



Ho = 65±9 kms -1 Mpc -1 



(6) 



which is in reasonable agreement both with the well-known 
WMAP result, 71±3.5 kms -1 Mpc -1 , and with the most 
recent Hu bble diagram for a complete sample of type la su- 
pernovae jReindl et al.l Eo05). whose interpretation requires 
a slightly reduced value of Ho. 

On the other hand, at the given ratio H { loc) /Ho we have 



Odo 
Oao 



Ho 



H, 



(loc) 



(7) 



So, using the most popular value of the total Hubble con- 
stant Ho = 71 kms - Mpc -1 and our value of the local Hub- 
ble constant //q 1oc ' =56 kms -1 Mpc -1 , we get Qdo/^ao ~ 
0.6, i.e. a much less fraction of the dark energy than the 
commonly-accepted one. Therefore, a slightly reduced value 
of Ho seems to be a preferable option. 



5 CONCLUSIONS 

As follows from the above analysis, the presence A-term can 
give us a reasonable explanation of the anomalous increase in 
the lunar orbit, consistent with the large-scale astronomical 
data. Thereby, this is one more argument in favor of the 
dark energy. 

Besides, if the local Hubble expansion really exists, it 
should result in profound consequences not only for cosmo- 
logical evolution but also for the dynamics of planetary sys- 
tems and other 'small-scale' astronomical phenomena, which 
have to be studied in more detail. 
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APPENDIX A: A POINT-LIKE MASS IN THE 
LAMBDA-DOMINATED UNIVERSE 

Solution of the equations of General Relativity for a point- 
like mass M in the Universe filled only with A-term is 



ds 2 



2GM Ar u 



+ 1- 



2GM Ar" 



c at 



dr 



+ r' 2 (d6 2 + sin 2 6 dip 2 ) (Al) 

llKottleil Il918f) . where G is the gravitational constant, 
and c is the speed of light (for general review, see also 
iKramer et"aTlll98fJ) . 

After a transformation to the cosmological Robertson- 
Walker coordinates 



ao exp 



t' 



t 



1 r 



hi 



1 



2 c 

metric HA It takes the form 



a 2 2ct\ 
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r V ro J 



ds 2 = gttcdt 2 + 2 g tr cdtdr + g rr dr 2 



+ gee d8 2 + g vv dtp 2 . 



where 
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(A6) 



Qrr — 



-1-^4- 
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(A7) 



(A8) 



gee = g^/sin 2 ^ = 
In the above formulas, r„ 



, „ - 2GM/c\ r = y/3/A , and ao is 
the scale factor of FRW Universe. 

Taking ao — 1 at t = and keeping only the lowest-order 
terms of r g and 1/ro, we get 



gw- 



gtrf 



2GM 



1 - 



VAi 

"7T 



agmVX 

V3c 2 



1 + 



2GM 



2c^fKt 

~7T 



gee = J w /sin 



(A9) 
(A10) 
(All) 
(A12) 



As is seen from the above expressions, manifestation 
of the A-term in some components of the metric tensor 
looks like the influence of variable G, if we assume that 
G = Go+Gt, where G — — cVA/-\/3- Unfortunately, such in- 
terpretation is not self-consistent: the A-dependence of other 
components is irreducible to the effect of variable G. 



APPENDIX B: STATISTICAL ANALYSIS OF 
THE EARTH'S ROTATION DECELERATION 

It is well known that diurnal period of the Earth Te, along 
with secular increase in the course of time, experiences 
a wide range of oscillations with (quasi-)periods from an- 
nual scale to many decades and, probably, even longer. In 
fact, just the long-period variations enforced the most of 
researchers to supplement the available telescopic data on 
the Earth's rotation by ancient records of the solar eclipses 
(whose reliability was, of course, not so good). 

On the other hand, as was already mentioned in the 
main body of the paper, we are interested in using only the 
data that are as close as possible to the LLR measurements. 
Of course, in deriving the secular term from the short data 
series, a special care should be paid to a probable inter- 
ference from the oscillatory components: if periods of such 
oscillations are comparable with the length of the series an- 
alyzed, they can partially mimic the secular (linear) term 
and, therefore, produce a substantial error in its final value. 
So, the basic idea of our analysis was to simulate such mimic 
effect and to find its maximum contribution to the secular 
term. 

Let the observational data be fitted by the following 
function: 



fit) = /o + fi t + a sin(27rf/T) + b cos(2nt/T) . 



(Bl) 



where T is the trial period; and fo, fi, a and 6 are the un- 
known regression coefficients. They are determined, as usual, 
by minimization of the functional 
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1656 1666 1676 1686 1696 N b 



Figure Bl. Values of (solid curve), j( mm ) (dashed curve) 

and j( max ) (dotted curve) as functions of the beginning of the 
data series N\>. In the intervals where some of the curves are not 
plotted, the corresponding extrema were not found. 



-We 

$(f ,fx,a,b,T,N h ,N e ) = 2 [/(**)- /if (B2) 

i=N b 

with respect to fo, fx, a and b. Here, i is the number of year, 
fi is the corresponding observed value of /; JVb and N e are 
the beginning and end of the analyzed time interval. 

The resulting values of fo, fx, a and 6 will be, in general, 
functions of T, N b and AT e - The extent of 'mimic' contribu- 
tion from the periodic components to the secular term fx 
(i.e. the contribution under the worst conditions) can be 
characterized by the quantities 

f[ min> (JV b , iV e ) = mm fx (T, N h , N e ) , (B3) 
/I™ 1 (N b ,N e )= max fi(T, N b ,N c ) ; (B4) 

while f[ hn \Nx,, N c ) will denote the value derived with the 
purely linear regression function 1B1L when the coefficients 
a and b were dropped out a priopi. 

The second important task is to estimate the effect of 
truncation of the time series on the resulting value of the 
secular term. (Since the earliest telescopic data might be 
not so reliable, it may be reasonable to exclude them from 
the analysis.) So, we performed a number of calculations 
with the shorter data series, when /Vb changed from 1 to 50 
(while N c was always the same). 

As was already mentioned, the primary o bservational 
data fj were taken from the monograph by ISidorenkovl 
(2002). They represent the mean annual variations in the 
Earth's angular velocity f(t) = (<5S1e/^e) X 10 10 for the 
years 1656 to 2000. 

The main results of our regression analysis are shown 
in Fig. IB1I As is seen, the secular term fx can vary in to- 
tal (as function of both T and A^b) from —1.135 yr _1 to 
—0.895 yr _1 . So, its average value, distorted by the 'mimic' 
effect of periodic variations and insufficient accuracy of the 
earliest data, can be written as fx = —1.015 ± 0.12 yr _1 . 
Then, the required rate of the Earth's diurnal deceleration 
will be T E = - T E 10" 10 /i (where T E = 8.64xl0 4 s), result- 
ing in T E = (8.77±1.04)xlO~ 6 syr" 1 . 



